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Abstract 

The production of pairs of hadrons in hadronic collisions is studied using a next-to-leading- 
order Monte Carlo program based on the phase space slicing technique. Up-to-date fragmentation 
functions based on fits to LEP data are employed, together with several versions of current parton 
distribution functions. Good agreement is found with data for the dihadron mass distribution. A 
comparison is also made with data for the dihadron angular distribution. The scale dependence 
of the predictions and the dependence on the choices made for the fragmentation and parton 
distribution functions are also presented. The good agreement between theory and experiment is 
contrasted to the case for single vr'^ production where significant deviations between theory and 
experiment have been observed. 
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I. INTRODUCTION 



The many successes of QCD at describing large momentum transfer processes have helped 
establish it as the theory of the strong interactions. Indeed, largely due to this success, 
research concerning QCD has moved from testing the theory to testing the approximations 
used to obtain predictions from the theory. Even though the overall description of large 
momentum transfer processes appears to be satisfactory, there are still some systematic 
discrepancies between the theory and experiment. These include, for example, problems 
observed in direct photon production [|T| and single tt^ production . One phenomenological 
approach has emphasized that the single particle production processes are sensitive to recoil 
corrections due to the emission of initial state radiation - also known as smearing 0. 
Another viewpoint [@] has stressed that, at least in the case of 7r° production, there may 
be problems with our knowledge of the fragmentation functions in the region where the 
momentum fraction, z, taken by the produced particle is large, since this is a region where 
the data to which these functions are fitted are limited. Also, it has been noted that this 
same high-z region may require significant threshold resummation corrections. 

The production of hadron pairs relies on the same underlying dynamics as single parti- 
cle production. Furthermore, the production of high-mass pairs relies on the same high-^ 
region of the fragmentation functions as does single particle production. Thus, if threshold 
corrections are important, or if the fragmentation functions are inadequately known, then 
one might expect to see comparable disagreement between the data for dihadron production 
and the predictions as one sees for single particle production. 

It is the purpose of this paper to present a comparison between data for high-mass 
dihadron production and predictions based on QCD. These predictions have been obtained 
using a next-to-leading-order Monte Carlo based program which uses a variant of the phase 
space slicing technique ^. This allows the same kinematic cuts used in the extraction 
of the data to be imposed on the theoretical predictions. Another treatment of dihadron 
production at next-to- leading-order using Monte Carlo techniques may be found in [§. In 
addition, the related process of photon pair production, including photon fragmentation 
functions, may be found in 0. The formalism in these two references has been used to 
investigate 7r°7 and tt^tt^ production at the LHC as reported in p]. 

The outline of the paper is as follows. In the next section a brief overview of the structure 
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of the calculation is presented. Then, the predictions of the program are compared to data 
in Sec. |T^. A summary and some conclusions are presented in Sec. 0. 

II. NEXT-TO-LEADING-ORDER MONTE CARLO CALCULATION 

The calculation described in this section is based on the two-cutoff phase space slicing 
technique describd in [||, |^. The basic concept is to partition the three-body phase space 
into three regions using two cutoff parameters, 6s and 6c- One region is where the 2-^3 
matrix elements have soft or coUinear singularities, one contains hard-collinear singularites, 
and in the remainder the matrix elements are finite. In the soft and hard-collinear regions 
the matrix elements are approximated using the soft or leading pole approximations, re- 
spectively, and the variables describing the soft or coUinear quanta can be integrated over 
analytically. The results have the same form as the lowest order 2—^2 contributions, 
but depend explicitly on the cutoffs used to define the the soft and hard-collinear regions. 
Likewise, the remaining finite 2^3 contributions depend on the cutoffs used to isolate 
the divergent regions. These two types of contributions are used to generate two-body and 
three-body weights which are added together at the histogramming stage. For infrared-safe 
observables the dependence on the cutoffs cancels, provided that sufficiently small values of 
the cutoffs are chosen. Specific examples of this procedure are given in [^. 

The case of high-mass dihadron production is formally rather similar to that for single 
hadron production with the addition of another fragmentation function. The treatment 
presented here, therefore, follows closely the presentation given in Sec. Ill E. of Ref. 0. 

The input needed for this calculation includes the squared matrix elements for the 2 — 3 
subprocesses and the results for the 0{al) one-loop contributions to the 2 — > 2 subprocesses 
0, |T0|. For the purpose of this example, the notation of [|T0| will be used, since much of 
the input needed can be found in the appendices of that paper. The partons are labelled 
as A + B ^1 + 2 and A + B ^1 + 2 + 3 for the 2 ^ 2 and 2-^3 subprocesses, 
respectively. A flavor label is used to denote the flavor of parton A, and similarly for the 
other partons. 

The lowest-order contribution to the inclusive cross section for producing two hadrons 
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hi and /i2 in a collision of liadrons of types A and B can be written as 
1 



da 



B 



X 



2xaXbs 
(47ra, 



GaA/A{xA)GaB/B{xB)Dhl/ai{zi)Dh2/a2{z2)dXAdXBdZi dz2 
aA,aB,ai,a2 



2.1) 



w{aA)w{aB) 

where a = {0^1,05,01,02} and p = {pa^Pb,PuP2} denote the sets of flavor indices and 
parton four- vectors, respectively. The factors appearing in the spin/color averaging are 
given by 

2{1 — e)V a=gluon 
2N a=quark or antiquark 

with = 3 and V = N"^ — 1. The factor dT2 is the differential two-body phase space element 
in n-dimensions, 

d^~^p^ d^~^Pi 



w{a) 



dT^ 



-(27r)"5"(pi+p2-P3-P4). 



(2.2) 



2j9||(27r)"-i 2pl{27iy 

Eq. ( p.l[ ) gives the contribution where parton 1 fragments into the hadron hi and parton 
2 fragments into /i2- Care must be taken to explicitly include in the sum over a those 
terms corresponding to the case where parton 2 fragments into hi and vice versa. For 
compactness, these terms will not be explicitly written. The squared matrix elements for 
the various subprocesses, denoted by ip^^^a^p), may be found in Ref. [|10|1 . 

Next, consider the one-loop virtual corrections to the 2^2 subprocesses. These take 
the form 

1 

2xaXbS aA,aB,ai,a2 



da'' 



GaA/A{xA)GaB/B{xB)Dh^/a^{zi)Dh2/a2{z2)dxAdXBdZi dz2 



where 



w{aA)w{aB) 



Til 



2tt \2pa-PbJ r(l-2e) 



tlj^^\d,p)dr2 



(2.3) 



(4)/ 



-^E^K)-7E7(c 



+ - y In 
2e ^ 



Pm ' Pn 

,Pa ■ Pb , 



- ^j:^^'\d,p)+^ljPsia,p) + Oie). 



(2.4) 



This expression for ip^^'^ differs slightly from Eq. (35) in Ref. |]T0| because we have chosen to 
extract a different e dependent overall factor: a factor of r(l -|- e)r(l — e) ~ 1 + has been 



absorbed into the above expression for ip^^\ Furthermore, the arbitrary scale Q^^s used in 
Ref. |]lOl has been chosen to be 2pA ■ Pb- The expressions for the functions V'm'n ^'^d ^/'^^ 



may be found in Appendix B of Ref. |0. The quantities C(a„) and 7(a„) are given by 

C{a) -- 



N = 3 a=gluon 
Cf = i a=quark or antiquark 



and 

f {UN -2nf)/6 a=gluon 
7(a) = <^ 

I 3Cf/2 a=quark or antiquark 

It will be convenient for subsequent expressions to adopt the following notation: 

■^-[2pa-pb) r(i-26)- ^'-^^ 

The one loop virtual contributions can now be written as 

0?Cr" = ^ GaA/A{xA)GaB/B{xB)Dh^/a^{zi)Dh^/a.,{z2)dxAdXBdZidZ2 

X ,Jl^^^f4i + iI + ^;).r, (2.6) 

where 



w{aA)w{aB) 2tc \ e"^ e 



Al = - Y.CM^^\d,p) (2.7) 

n 

A\ = -E7(an)^(^n«,Pl + ^ E In f^^l i^l^fr^{S,P} (2.8) 

V ^ \pa-pbJ 

^0 = -^Y.C{an)^P^'\3,p) + ^pPs{S,P). (2.9) 



Next, the contributions from the 2 — 3 subprocesses in the limit where one of the final 
state gluons becomes soft are needed. The contributions of the 2 — > 3 subprocesses may be 
written as 

d(^^^^ = TT— /^T'^^ \^ GaA/A{xA)Ga/B/B{xB)Dh,/a^izi)Dh^/a:,{z2) 

2XAXBSW{aA)W{aB) aA,aB,a^,a2,az 

X '^{aA,aB,ai,a2,a3,PA,PB,PuP2,P3) dT3dxAdxBdzidz2. (2.10) 

The expressions for the 2-^3 squared matrix elements appearing in Eq. ( |2.10|) may be 
found in Ref. M. As noted earlier for the two-body contributions, one must include in the 
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sum all possible parton to hadron fragmentations. Here dT^ is the three-body invariant 
phase space factor in ra— dimensions: 

d'^'^p^ d^~^Pi d'^~^p5 



dr. 



^p^(2vr)"5"(pi +P2-PS-P.- p,) . (2.11) 



2p||(27r)"-i 2p^(27r)"-i 2pl 
Consider the case where the soft gluon is parton 3. In this limit, the function ^ may be 
expanded as: 



Pm ' Pn 



m<in 



Pm ■ PsPn ■ P3 



-^Inln (flA , flfi , ai , ^2 , Pa , , Pi , ; 

(2.12) 



Next, one must integrate over the soft region of phase space defined by < 
6s\J1pA ' Pb/'^- This is easily done using the integrals given in the appendix of Ref. [^. 
The resulting soft contribution may be written as 
1 



da' 



2xaXbs 
(47ra. 



XI GaA/AixA)Gas/BixB)Dh^/a^ {zi)Dh.,/a2 {z2)dXA dxs dzi dz2 
aA,aB,ai,a2 

■J'^(^ + ^ + A'o]dT2 (2.13) 



w{aA)w{aB) 2tt \ e 



where 



= J2C{a,,)ij^'\d,p) 

n 

Al = -21n5,5:C(o„)^(^)(o,p-) - i ^ In f ^ 

n ^ m,n \PA 

= 21n2 5,XC(o„)V^^'^(a,p) 



Pb 



(2.14) 
(2.15) 



+ 



+ 



(4,c) , (4,c) 
A,l + rB,2 



\Pa-PbJ \Pa-Pb, 



+ 2 In (5^ In 



Pa-Pb, 

P2-P3^ 



.(2.16) 



\Pa ■ Pb J \Pa ■ Pb J \Pa ■ Pb , 

After the coUinear singularities associated with the parton distribution functions and frag- 
mentation function have been factorized and absorbed into the corresponding bare functions, 
there will be soft-collinear terms left over due to the mismatch between the integration lim- 
its of the collinear singularity terms and the factorization counterterms. Collecting together 
the various collinear terms, the result can be written as follows: 
1 



da 



coll 



X 



2xaXbs 
{Anas 



w{aA)w{aB) 27r 



GaA/A{xA)GaB/B{xB)Dh^/aAzi)Dh^ia,^{z2)dXAdXBdZi dz2 

(2.17) 



aA,a.B,ai,a2 
\2 ^, / A coll 



+ Af' ) dT2 



where 



Af' = Y.[2\n5sC{an)+l{an)] (2.18) 

n 

= [21n5.C(a„)+7(a„)]lnf^^^') 

n=A,B \ / 

+ ^j21n5.C(aO+7K)]ln(^^^). (2.19) 

Here /i/ and Mj are the initial and final state factorization scales. 

After the mass factorization has been performed, the bare parton distribution functions 
and fragmentation functions have been replaced by scale dependent MS functions. In ad- 
dition, there are finite remainders involving functions G and expressions for which may 
be found in 

da = — 7^ — — r—il)^'^\d,f))dxAdxBdzidz2dT2 

2xaXbS aA,a^a^A2 W[aA)w{aB) 2tT 

+ GaA/A{xA, IJ,'j)GaB/B{xB, IJ^})Dh,/aj (Zi, Mf)Dh^/a^ {Z2, Mf) 
+ Ga^/AixA, l^})GaB/BixB, IJ,})Dh,/a, (^1, M'j)Dh^la^ (^2, Mj) 

+ Gaji/AixA, IJ^'j)GaB/B {xB,^i})Dh,/aM,Mf)Dh,/a,iz2,Mf)\ . (2.20) 

At this point, all of the singular terms have been isolated as poles in e or have been 
factorized and absorbed into the bare parton distribution and fragmentation functions. The 
e dependent pole terms all cancel amongst each other: 

A^ + A^ = (2.21) 
Al + Al + Af^ = 0. (2.22) 

The finite two-body contribution is given by 
rfcr^-"2 = da^ + da 

+ —r—^l-^^^A/AixA, fJ'})Gas/BixB, l^}) 

XaXbS aA,a^a^,a2 w{aA)w{aB) 

X Df,^/a, (^1, MhDh,/a, (^2, M?) ^ \AI + AI + Af'] dxA dxB dzi dz2 dT2.{2.23) 

The three-body contribution, now evaluated in four dimensions, was given in Eq. (p. 10 ) 
where now the soft and coUinear regions of phase space are excluded. 
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The structure of the final result is two finite contributions, both of which depend on the 
soft and collinear cutoffs - one explicitly and one through the boundaries imposed on the 
three-body phase space. However, when both contributions are added while calculating an 
observable quantity, all dependence on the cutoffs cancels when sufficiently small values of 
the cutoffs are used. 



III. COMPARISON TO DATA 

Two sets of next-to-leading-order fragmentation functions have become available recently 
||TT| , Both sets have been fit to high statistics data from e^e~ experiments. Accordingly, 
only charge symmetric combinations, e.g., + h~ , have been determined and the sets do 
not have fragmentation functions for individual charge states. Nevertheless, these sets can 
be used to generate predictions for experiments which measured either vr^vr^ final states or 
symmetric combinations of charged hadrons. The NA-24 []13[, CCOR iH], and E-706 |T5 



experiments each measured the production of 7r° pairs while the E-711 ||I6[ experiment mea- 
sured the production of h~^h^,h~h~, and/i^/i^ pairs. In the latter case, one can combine 
the tabulated results to give the cross section for producing the symmetric combination 
{h+ + h-) + {h+ + h-). 

In the following, unless otherwise stated, the theoretical results have been obtained using 
the CTEQ5M [jl^ parton distributions and the KKP fragmentation functions. For the 



calculation of the cross section at fixed values of the dihadron mass,M, the renormalization 
and factorization scales have been chosen to be proportional to M, as this is the only 
observed hadronic variable with the appropriate dimension. 

In FIG. |1| the E-711 ||I6| data are shown with the leading-order (LO) and next-to-leading- 
order (NLO) theoretical results with two choices for the common factorization and renormal- 
ization scales. Cuts were applied to the rapidity of the pair, Y, the transverse momentum 
of the pair, pxpair, and cos^*, an estimate of the cosine of the scattering angle in the parton- 
parton center of mass frame. For these data, cos 6* was defined by first transforming to the 
frame where the momentum of the hadron pair had no longitudinal component. In general, 
the two hadrons will not be exactly back-to-back in this frame, due to their differing values 
of transverse momentum. The two values of the cosine of the angle between the hadron 
direction and the beam direction were averaged to obtain cos^*. The cuts used for the data 
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FIG. 1: Comparison of the NLO (solid) and LO (dashed) results with data from the E-711 experi- 
ment [p^]. 

shown in FIG. |I] were -0.4 < Y < 0.2, pTpair < 2 GeV and | cos^*| < 0.25. The NLO 
results can be seen to bracket the data while, for the scale choices shown, the LO results 
are significantly below the data. The large scale dependence evident at lowest order is due 
to the two powers of in combination with the scale dependence of the four distribution 
and fragmentation functions. In the kinematic regime covered by the data, the distribution 
and fragmentation function momentum fractions are large, so that the functions decrease 
with increasing values of the scale. These scale dependences result in a significant decrease 
of the cross section with increasing scale, as shown in FIG. |1]. The band covered by the 
corresponding NLO curves is narrower, although significant scale dependence remains. This 
is further examined in FIG. where the NLO results are compared to the E-711 data for 
three choices of scale. The format (data-theory) /theory is used in order to more clearly 
show the scale dependence. 

In FIG. I the depend ence of the results on variations in the choice of parton distributions is 
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FIG. 2: CTEQ5M results compared to the E-711 data for different scale choices. The renormal- 
ization and initial and final state factorization scales have been set equal to each other. 
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FIG. 3: Comparison between the CTEQ5M predictions with those from the CTEQ5HJ |17|set and 
several MRST H] sets. 
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FIG. 4: Relative contributions of the qq,qg, and gg subprocesses to the CTEQ5M results with 
H = 0.85M. 

shown, relative to the CTEQ5M results. For each curve the scale has been set to fi = 0.85M. 
The CTEQ5HJ set has a gluon distribution which has been enhanced in the high-x region 
in order to better describe the high-i^r jet data from the CDF and D0 Collaborations. 
Relative to the CTEQ5M distribution results, one can see an overall increase in the cross 
section, with the increase becoming larger towards the high mass end. Note that an increase 
in the scale from 0.85M to M would bring the CTEQ5HJ curve down to the level of the 
data. Also shown are the results for two of the MRST sets, one with the standard gluon 



and one with a reduced gluon. Although the two curves lie below the CTEQ5M results, a 
modest decrease in the scale choice would raise the curves to be in accord with the data. The 
results for a third MRST set with an increased gluon distribution are essentially identical 
with the CTEQ5M results. 

In FIG. ^ the relative contributions of the quark-quark, quark-gluon, and gluon-gluon 
subprocesses are shown versus ^/T = M/ ^/s. This dimensionless variable is approximately 
the value of the parton momentum fraction which is probed in the production of the high 
mass hadron pair. As expected, quark-quark scattering dominates at the upper end of the 
mass range covered by the data. Nevertheless, there is a significant contribution from quark- 
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FIG. 5: Comparison between the results obtained using the KKP [11| and BFGW |12] fragmenta- 
tion functions. 

gluon scattering over much of the mass range. This is similar to the situation for high-ii^y jet 
production and, indeed, the results in FIG. ^ do show some sensitivity to the choice of the 
parton distributions, e.g., CTEQ5M vs. CTEQSHJ. Unfortunately, the scale dependence, 
even at NLO, is such as to preclude favoring one set over the other. 

Next, in FIG. |^ the dependence on the choice of the fragmentation functions is shown. 
The results from the two sets agree to within about 10% across the mass range shown. 

From the results shown thus far, several conclusions can be drawn. First, the NLO results 
give a very good description of the mass distribution for symmetric hadron pairs measured by 
the E-711 Collaboration. The variations observed due to different choices of the distribution 
and fragmentation functions are easily compensated for by changes in the renormalization 
and factorization scales. Nevertheless, extreme variations of these scales are not needed in 
order to describe the data. 

Next, consider the data for producing pairs of neutral pions. In FIG. ^ the theoretical 
results are compared to data for the process pp n^-n-^ + X as measured by the NA-24 |l^ 
and CCOR 0] Collaborations. The cuts used for these data are —0.35 <Y< 0.35, pxpair < 
1, and I cos^*| < 0.4. The same scale choice of fi = 0.85M as used for the E-711 data gives 



12 



X, 
PP — > 71 71 + X 

CTEQ5M n = .85 M 



10 



CCOR Vs = 62.4GeV 
CCOR Vs = 44.8GeV 



> 



> 

T3 

D 
T3 



0.1 



0.01 . 



0.001 . 



0.0001 




M (GeV) 

FIG. 6: Comparison between the CTEQ5M results and data from the NA-24 [||] and CCOR 
CoUaborations. 

a good description of the experimental results. Again, the CTEQ5M and KKP distribution 
and fragmentation functions have been used. The same results are shown in FIG. in a 
data/theory format versus ^/T. For the NA-24 data, the statistical and systematic errors 
have been added in quadrature. No discussion of errors was contained in the CCOR paper 
\L^ . However, in an earlier publication on single tt^ production the CCOR Collaboration 



quotes an overall error of 25% for both energies with an additional 5% relative normalization 
error between the results for the two energies. Thus, the three data sets are seen to agree 
within the quoted errors. Furthermore, given these errors, the agreement with the theoretical 
results is acceptable, although deviations are apparent in FIGS. ^ and^. 

Noting the tendency for the data to lie below the theory for small i/r and above it 
for larger ^/T, it is natural to ask whether the CTEQ5HJ distributions might yield better 
agreement with the data. These results are shown in FIG. ^ with a scale choice of fi = M. 
Although the agreement is slightly better, the modified gluon distribution is unable to bring 
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FIG. 7: The same as for FIG. |^ presented in the data/theory format. 
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FIG. 8: The same as for FIG. except using the CTEQ5HJ distributions. 

the theory and the data into complete agreement. Nevertheless, given the size of the errors, 
no definitive conclusion can be drawn. 
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FIG. 9: Comparison between the NLO results and data from the E-706 Collaboration |15|. 



In FIG. ^the NLO predictions are compared with data for dipion production as measured 
by the E-706 Collaboration |T^. For these data the pions were separately required to satisfy 
Pt > 2.5 GeV and -0.8 < y < 0.8 (-1.05 < y < 0.55) for the ^ = 31.6 (38.8) GeV 
data. The difference between the azimuthal angles for the two pions, A0, was required to 
be greater than 105 degrees. No cuts were placed on cos 6** or on pxpair- There is good 
agreement between the theory and the data when a scale choice of fi = 0.35M is used. This 
value is significantly smaller than that used in the previous comparisons. Of course, the 
cuts used for the E-706 data set are far different than those used for the other sets. As a 
consistency check on their data the E-706 Collaboration [|1^ also presented results obtained 
using cuts similar to those used for the E-711 data. These data are compared to the NLO 



predictions in FIG. |10| with fi = 0.50 M and 0.85 M. The two curves bracket the data, 
suggesting that E706 results are compatible with those of the other experiments when the 
same cuts are used. This suggests that the need for a much smaller scale when comparing 
with the E-706 data in FIG. is a problem due to the calculation not being able to properly 
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FIG. 10: Comparison between the NLO results and data from the E-706 Collaboration [^]. 

reproduce the effects of the different sets of cuts. 

In order to investigate this situation further, consider the effects of a cut in A0. The 
signifigance of the A0 cuts lies in the fact that the observed A0 distributions are broader 



than those given by the NLO calculations, as can be seen in FIGS. |ri| and |T2|. At lowest 
order, with coUinear fragmentation and distribution functions, the theoretical predictions 
have the two hadrons being produced back-to-back, i.e., with A0 = 180 degrees. At next- 
to-leading-order, the 2 — i> 3 subprocesses allow for the A0 distribution to develop a non-zero 
width. Nevertheless, it is still narrower than the experimental observations. The acceptance 
for the CCOR and E-711 experiments was limited to A0 > 140 degrees, and this cut was 
also placed on the NA-24 data shown previously, so essentially none of the NLO generated 
cross section was rejected. Opening up the cuts to the value A0 > 105 degrees used by the 
E-706 experiment does not, therefore, result in any increase of the theoretical cross section. 
However, it does result in an increase in the experimental cross section. The net result is 
that the NLO calculation will not correctly reproduce the effects of different Acp cuts because 
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FIG. 11: Comparison of the NLO A4> distribution with data from E-706 [15| 
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FIG. 12: Comparison of the NLO A(j) distribution with data from E-706 Jll. 



the theoretical distribution is too narrow. This results in a relative normalization shift when 
comparing experiments which used different cuts in A0. This shift can be accomodated 
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by altering the renormalization and factorization scales used in the theoretical calculations. 
Alternatively, one can restrict the comparison to data sets which use the same A0 cut. 
Similar considerations apply to cuts on other variables such as the net transverse momentum 
of the hadron pair, pxpair- The comparison with E-706 data is shown in FIG. where 
one can see that the theoretical pTpair distribution is significantly narrower than is observed 
in the data. Comparisons to data which integrate over the full pxpair distribution will differ 
from comparisons to data sets which place cuts on this variable. 

This situation should, in fact, come as no surprise. Both the A0 and prpair distributions 
are examples which are delta functions at lowest order. Non-trivial contributions to these 
observables only start in next-to-leading order. In that sense, the curves shown here for these 
distributions are really leading-order only and they diverge at the endpoints corresponding 
to the third parton being soft and/or coUinear (A0 = 180° or pxpair = 0). A more realistic 
treatment of these observables would require the application of soft gluon resummation 
techniques. Note, however, that there are compensating singularities at the endpoints of the 
distributions, so that a finite result is obtained after integrating over the distribution and 
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FIG. 14: Comparison between the CTEQ5M results and the E-711 angular distribution data for 
7 < M < 7.5 GeV using two choices for the scale parameter. 

the normalization of the integrated distribution is thus calculated to next-to- leading order. ^ 
From the standpoint of comparing to NLO calculations, it would be better if the data 
sets were defined only by cuts on variables whose distributions are well described by the 
calculation. In this sense, the E-706 procedure is to be preferred since a large portion of 
both the A(j) and prpair distributions were integrated over. 

The cuts utilized in the analysis of the E-706 data differ substantially from those which 
were used for the CCOR, NA-24, and E-711 data sets. Therefore, separate comparisons are 
required and the optimum choice of scale will differ between the two sets of experiments. 
It must be stressed that this is not an experimental problem, but rather the result of the 
fact that the NLO calculation does not properly describe the distributions in some of the 

^ Due to the use of finite width bins, the divergent behavior at the endpoint of the pTpair distribution in 
FIG. O is not evident. The first bin remains finite as it contains the endpoint contribution, as well. 
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variables used for making the cuts. 

The E-711, E-706, and CCOR Collaborations each measured the angular distribution of 
the dihadron pair in the parton-parton center-of-momentum frame. None of the experiments 
observed any significant variation of this distribution with dihadron mass or with energy. 
The theoretical results for the normalized cos 6** distribution are compared to the E-711 data 



[|T|, H in FIG. [T| for 7 < M < 7.5 GeV and in FIG. |T| for M > 7.5 GeV. In each figure 
two curves are shown corresponding to the NLO results with scale choices of 0.85M and 



< Pj^ >, where < p"^ > is the average of the squared transverse momentum for the two 
observed hadrons in the event. Note that for the case of two-body kinematics at fixed M, 
the parton transverse momentum is ^ sin 9* and the parton and hadron transverse momenta 
are nearly the same since the fragmentation variable z is near one. Thus, one can argue 



that for fixed M and 9* either M or J < pj^ > is a. valid choice for the scale. The choice of 



'< > gives a steeper distribution which is in better agreement with the data than is the 
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result obtained with the choice of 0.85M for the scales. This steepening occurs because at 
fixed M as cos^^* 1, ^J < p\ > decreases. This decreasing scale causes an increase in the 
theoretical cross section. However, since the distribution is normalized to unity at cos^* = 0, 
the end result is a steeper angular distribution. 

IV. SUMMARY AND CONCLUSIONS 

A next-to-leading-log Monte Carlo program has been constructed for symmetric dihadron 
production using the two cut-off phase space slicing formalism described in Refs. [|, Q . This 
process serves as a probe of the underlying hard-scattering subprocesses which complements 
that provided by single particle production. For high mass pairs the relevant range of the 
fragmentation momentum fraction z is comparable to that for single particle production. 
The results presented here show that the NLO QCD formalism is capable of giving a good 
description of the data for the dihadron mass and angular distributions. There appears 
to be no anomalous behavior with respect to either the dihadron mass or the center-of- 
mass energy. This is in marked contrast to the cases for direct photon and single hadron 
production at fixed target energies. This process, therefore, provides encouraging evidence 
that the underlying hard scattering is correctly described by QCD and that the problems 
with the single photon and single hadron cross sections may be ascribed to a combination of 
of effects due to an incomplete application of the theory and possible inconsistencies amongst 
the various data sets. 
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